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ABSTRACT 

We detect a significant broadening in the wings of the 401 Hz peak in the 
power spectrum of the accreting millisecond binary pulsar SAX J1808. 4-3658. 
This feature is consistent with the convolution of the red noise present in the 
power spectrum with the harmonic line. We conclude that the flux modulated 
by the spin period shows aperiodic variability similar to the red noise in the 
overall flux, suggesting such variability also originates at the magnetic caps close 
to the neutron star surface. This is analogous to the results found in some longer 
period, higher magnetic field pulsators in high mass X-ray binaries. 



Subject headings: pulsars: individual (SAX J1808. 4-3658) — accretion — stars: 
magnetic fields — stars: neutron — X-rays: stars 



- 2 - 



1. Introduction 

The power density spectra of low mass X-ray binaries show broad band noise com- 
ponents, produced by aperiodic variability, the sum of which typically rises towards lower 
frequencies (sometimes flattening below a frequency z/ knee ) and is often termed red noise. 
Other power spectrum features are often present, including relatively broad peaks, quasi 
periodic oscillations, bumps, wiggles and low frequencies excesses (see e.g. Wijnands & van 
der Klis 1999 for a review). In many respects these noise components are similar to those 
observed in high mass X-ray binaries, containing X-ray pulsators and black hole candidates 
(see e.g. Belloni & Hasinger 1990). The continuum component of the red noise is usually 
modeled by broken power laws and/or Lorentzians. 

In power spectra of X-ray pulsators a number of harmonic lines arising from the mod- 
ulation of the flux due to the neutron star spin is also present. In this case, information 
on the regions giving rise to the aperiodic variability can be extracted from the study of 
the power spectra (Burderi et al, 1997, hereafter B97; Lazzati and Stella, 1997, hereafter 
LS97). Indeed, if all or part of the aperiodic variability is produced by matter tunneled by 
the magnetic field lines and accreting onto the magnetic polar caps of the neutron star, a 
modulation of that part of the aperiodic signal at the spin period would be expected (har- 
monic coupling). In the power spectra this would appear as a broadening of the wings of 
the harmonic lines due the convolution of each harmonic line with the red noise (Burderi, 
Robba & Cusumano 1993) that produces a rescaled version of the red noise placed at both 
sides of each harmonic peak. Alternatively, if none of the aperiodic variability is generated 
in the tunneled matter, the periodic and red noise features are expected to be independent 
of each other and there should be no broadening at the base of the harmonic lines. 

A significant broadening in the wings of the harmonic lines was detected in the power 
spectra of some X-ray pulsators (B97; LS97), indicating that aperiodic variability is modu- 
lated at the spin period of the neutron star, i.e. produced at the neutron star magnetic polar 
caps. More recently, a coupling between the period and aperiodic variability of Hercules X-l 
has been reported (Moon & Eikenberry 2001). 

In 1998 the first millisecond accretion powered X-ray pulsar was discovered with the 
Rossi X-ray Timing Explorer (RXTE) by Wijnands & van der Klis (1998a). The pulsation 
frequency was ~ 401 Hz and showed a periodic modulation at 2.01 hours testifying to the 
presence of a low mass companion orbiting the neutron star with a mass function of the 
secondary star of 3.8 x 1O~ 5 M (Chakrabarty & Morgan 1998). The source was identified 
with SAX J1808. 4-3658, a transient first detected with Beppo SAX in 1996 which displayed 
type I X-ray bursts and was classified as a low mass X-ray binary at an estimated distance 
of 4 kpc (in 't Zand et al. 1998). This distance has been recently updated to 2.5 kpc (in 't 
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Zand et al. 2001). The coherent pulsation was detected during the 1998 outburst when the 
source was repeatedly observed with RXTE from the 11 th of April to the 6 th of May. The 
source reached a peak luminosity of 3 x 10 36 erg s _1 on the 13 th of April and then faded 
gradually below 10 35 erg s _1 by the Q th of May. Throughout this period the X-ray spectrum 
remained stable, featuring a power law shape with a photon index ~ 2 in the 15-100 keV 
band and a high energy cut-off above about 100 keV, similar to the X-ray spectra of other 
type I bursters in their low state (Gilfanov et al. 1998; Heindl & Smith 1998). Wijnands 
& van der Klis (1998b) studied the rapid aperiodic variability and revealed the presence of 
a red noise component in the power spectra flattening below a knee frequecy i^knee ~ 

1 Hz. 

The periodic modulation was nearly sinusoidal in shape and the power spectra displayed a 
very pronounced coherent peak at ~ 401 Hz. 

The possibility of detecting a modulation at the spin frequency in the aperiodic vari- 
ability depends crucially on the spin frequency: if this is too low, i.e. u sp - in <C u^ ncc the 
broadening of the harmonic line's wings is so wide that it blends in the red noise itself (B97; 
LS97). However, for SAX J1808. 4-3658 u sv - in ^> u^ ncc , which means that the typical duration 
of the aperiodic variations is many spin cycles, making this source a particularly well suited 
candidate to search for the presence of a coupling between this red noise component and the 
periodic variability. 

2. The model 

Analytic derivations of the power spectrum of a signal consisting of aperiodic variability, 
a periodic modulation and a coupling of the two was described by B97 and LS97. Our aim is 
to use an analytic function to fit the power spectrum of the SAX J1808. 4-3658 to derive the 
entity (if present) of the coupling. In line with B97, we adopt randomly occurring shots as 
a convenient description of the aperiodic variability giving rise to the red noise (shot noise 
model). In this model the aperiodic variability is produced by a random superposition (in 
the light curve) of shots. These instabilities might originate in the accretion flow onto the 
compact object. We note that we can apply this model without loss of generality, as the 
coupling effect is independent of the mathematical model adopted to describe the red noise. 
In the following we briefly summarize some useful definitions and results from B97. 

Our shot noise process Mga.it) is obtained from the random superposition of identical 
shots occurring at a constant mean rate A: the shot magnitude is S — J h(t)dt, where h(t) is 
the intensity profile of a shot. This means that the average intensity of the shot component 
is Tgn = XS. We indicate the power spectrum of the individual shot as \H{u)\ 2 . 
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A signal I(t), as seen by a distant observer, emitted by a rotating neutron star of spin 
period P is, in general, composed by: 

a) a (constant) background component Jbck; 

b) an emission (diffuse) arising from matter not funneled on the neutron star magnetic caps 
Jdf(£) and therefore not affected by the periodic modulation induced by the neutron star 
rotation; 

c) a modulated (localized) component arising from matter accreting onto the magnetic poles 
Ihc(t)- The periodic modulation function M(t), representing the geometric lighthouse effect 
produced by the neutron star rotation, is a series of identical pulse profiles of period P = 1/vo- 
Expanded in Fourier series M(t) = C + J2k c k cos(27r k v t + (f) k ). In the case of SAX J1808.4- 
3658 the pulse profile is almost sinusoidal and we can approximate M(t) ~ C + c\ cos(27r vot). 
We require < M(t) < 1, this poses a constraint on C and ci, namely, C = 1 — c±. Thus we 
have 

/(*) = 4ck + W*) + he (t) x M(t). 

Let us assume that the diffuse and localized emission are each composed of a shot noise 
component and a constant component: we assume that the shot magnitude S and profile 
h(t) is the same for both the shots arising in the localized and in the diffuse regions, while 
the shot rates are respectively Alc and Adf- Kdf and K LC are respectively the constant 
emission of the diffuse and localized components of the signal. Thus we have: 

I BF (t) = K DF + Af SQ (t) DF 

hc(t) = K LC + Kn(t)LC 

and for the mean total intensity: 

I = J bck + (K DF + X DF S) + (K LC + A LC S) x (I - ci) 

The power spectra of the signal described above, taking into account the finite duration 
T of the light curves and adopting the normalization of Leahy (1983) is: 

P(u) =P + P RN + P H L + P C PL + PWN (1) 

The first term is the power at zero frequency and is only dependent on the total intensity of 
the signal. 

P = NI 2 TW^(u) (2) 

where = 2/7 is the normalization constant according to Leahy and where W^(u) = 
sin 2 (7rz/r)/(7rz/r) 2 is the window function. We note that P = W fc = k/T where the 
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Fourier frequencies are computed at k — 1, (T/2Ai) — 1. For this reason P has not been 
included in the fit discussed in §3. The red noise (RN) term is 



RN 



N x (A DF + (1 - Ci ) 2 A L c)5 2 /rn(^) (3) 



where, in the simplest case of a shot noise process obtained from the random superposition 
of identical shots, /rn(^) = \H(v)\ 2 /MXX.{\H(i>)\ 2 }, in a more general case /rn(^) is a non 
negative function normalized to 1 that describes the shape of the red noise. The harmonic 
line term is given by: 

Phl = N x ( Cl /2) 2 x (K LC + X LC S) 2 x TW 2 (u - u ) (4) 

where, because of the almost sinusoidal shape of the SAX J1808. 4-3658 pulse profile, we have 
taken into account only the fundamental frequence. The coupling term is: 

Pcpl = N x ( Cl /2) 2 x \ LC S 2 f RN (v - u ) (5) 

It consists in a rescaled versions of the red noise placed at the harmonic frequencies, the 

scaling factor is r = Pcpl/Prn = \ ( j^r^ j + (i-c^Alc ) ' ^ ne wn ^ e n °i se component 
Pwn, induced by counting statistics, has, adopting the normalization of Leahy, an expected 
value of 2. The actual value is usually (slightly) less than 2 because of detector dead-time 
effects (e. g. van der Klis 1989), thus its value will be determined by the fit of the power 
spectra, as described in §3. 

We define the coupling fraction 77 = Alc/(Alc + Adf) as the ratio of the arrival rates 
of the localized shots over the total shot rate. Furthermore we define the pulsed fraction 
V f = (M(t) M Ax - M(t) M m)/M(t) M Ax = 2ci, where M(t) M Ax and M(t) M m are respectively 
the maximum and the minumum of the modulating function over one period. With these 
definitions < V f < 1 and we have ci/(l - Ci) = V f /(2 - Vf). Moreover, as M(t) > 
always, we have 1/2 < (1 — ci) < 1 and 1 < 1/(1 — ci) 2 < 4. Therefore the following relation 
must hold: 

" —2 ^ V f ^ 7 T ,-,-1/2 ( 6 ) 



1+ - , 



l-T) 



2v^ 1 + 4 (i=2) +2^ 



In the following section we will determine r from the fit of the power spectra and with 
relation (6) we will constrain the values of the pulsed fraction in relation to the fraction of 
coupled shots. 

Moreover we will be able to derive a relation between the diffuse component of the 
signal, as defined above, and the pulsed fraction. Averaging the signal only over the shot 
noise process we define: 

(/(*)) = /bck + (^df + XdfS) + (K hC + XlcS) x (C + ci cos(2tt u t)) (7) 
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thus 

(I{t))uAx = 4ck + (^df + A DF ,S) + {K hC + XlcS) 
(I(t))um = 4ck + (^df + XdfS) + (K LC + A LC .S) x (1 - 2 Cl ) 
We can therefore define the "observed" pulsed fraction 

V fohs = «/(*))max - (/(*))min)/((/(*))max " 4ck) 

with a little algebra we can derive V fohs = 2ci/ ( jggjaggf + l) = 2ci/ + l). Recol- 

lecting that = 2ci we can derive #4 = - 1. We can now define the diffuse fraction 
as the ratio of the diffuse component over the sum of the diffuse and localized components 
as f = i , — r. Using the relations above we have f = 1 — Vf ~ ha . We will see in the next 

^ (Jdf) + (-<LC) ° H 

section how the constraints on the pulsed fraction will have a bearing also on the fraction of 
the diffuse component. 



3. Data analysis and results 

We used the observations of SAX Jf808. 4-3658 performed by the Rossi X-Ray Timing 
Explorer (RXTE) during the April-May f998 outburst (see Table f). In particular we ex- 
amined the data from the Proportional Counter Array (PCA) (Jahoda et al. 1996) which is 
composed of a set of five xenon proportional counters sensitive to X-rays in the 2-60 keV 
range with a total area of ~ 7000cm 2 and a 1° FWHM field of view. We used data collected 
with ~ 122/is time resolution and 64 PHA channels which were available for all observations 
(except for the April 13 observation which was excluded from our analysis). 

We first corrected photons arrival times for satellite motion with respect to the solar 
system barycenter using the BeppoSax position for the source (in 't Zand et al. 1998) 
R.A. = 18/i08m29s, Dec = -36° 58.6' ( J2000). We then corrected for the orbital motion of 
the neutron star by using the orbital parameters derived by Chakrabarty & Morgan (1998). 
For our analysis we rebinned the data into bins of At = 488.28125;us (corresponding to a 
Nyquist frequency of 1024 Hz) and we used the events from all 64 PHA channels. The 
resulting lightcurves were used for the subsequent analysis. 

For each observation, we calculated consecutive FFTs over time intervals of T = 64s 
corresponding to power spectra with a frequency range of 0.015625-1024 Hz. These power 
spectra were then averaged to produce a single power spectra for each observation, 16 in all. 

The source flux decreased across the observations by a factor of ~ 100 (Gilfanov et 
al. 1998 ). The aperiodic variability was determined to change with flux level (Wijnands 
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& van der Klis, 1998b). On the other hand, the simple model of the aperiodic variability 
described in the previous section assumes that the characteristics of the signal do not vary 
with time. Thus in our analysis we selected and analyzed only the observations in which the 
source count rate was highest (between ~ 800 and ~ 500cts/s). In particular we selected 
the 6 consecutive observations between April 11 and 20. Examining the power spectra of 
each observation we consistently find an almost constant value for the red noise intensity 
and shape. 

The shape of the red noise in SAX J1808. 4-3658 was found to be reasonably well de- 
scribed by the sum of three empirical functions: 

3 



where: fi{v) = [1 + (27rTjZ/) ai ] 1 and pi indicates the normalization of each of them. 

In order to fit the spectrum efficiently, we logarithmically rebinned the spectrum so 
that the number of original bins averaged in the i th new bin follows n-i = int(1.02 l ), where 
int(x) is the integer part of x. The region around the coherent modulation peak (400.8 < 
v < 401.2) was left at the original Fourier resolution, while the intervals 352.37 < v < 400.8 
and 401.2 < v < 449.64, in which the peak broadening would be more important, were 
(uniformly) rebinned by a factor of 31 (see Fig. 1). The region around the second harmonic 
of the modulation was excluded (see above). We used a fitting function : 



where g(u) = W^(u) and f(u) = Pif\{i')+p 2 f 2 {i')+p 3 f 3 {i') as defined above. Thus Aif(v) = 
Prn, A 2 g{v - v Q ) = P HL , A 3 f(u - u ) = P CPL and A A = P WN . 

In order to check the robusteness of our fit we compared different fitting procedures, a) 
In the first place we have tried a global fit, that is we simultaneously fit the whole range 
leaving all the parameters of the fit free to vary. Then, keeping fixed the parameters that 
describe the three empirical functions f\, namely Pi, ati, Ti, we compare the fit on the whole 
range leaving A\, A 2 and A 3 free and leaving just A\ and A 2 free and forcing A% = 0, 
that is imposing the hypothesis of null coupling. We obtain two values of the x 2 with 
Ax 2 = Xi ~ x\ = 19- b) In second instance we repeated the procedure we compare the 
fit before and after freezing A 3 = only in the region near the armonic line, beetween 396 
and 406 Hz, the difference in the x 2 we obtain is A% 2 = 21. c) Last we wanted to check if 
the effect was still present freezing also A 2 , measuring the height of the harmonic line, thus 
we compared the fit in the same region of case b) and found that the A% 2 = 21. We thus 
conclude that the region that affects the fit is only the narrow region around 401 Hz and a 




i=i 



F(u) = AJ(u) + A 2 g(u - u ) + A 3 f(u - u ) + A A 
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significant coupling is present independently of the detailed fitting procedure. The values 
determined for case b) of the f\ and Ai, A 2 , A 3 parameters are listed in Tab. 2 along with 
their errors. From the best fit parameters A 3 and A\ we derive r = Pcpl/-Prn — 0.005. 
An F-test calculates the probability that the A% 2 variation is due to a chance improvement 
P = 8.38 x 1(T 6 . 

With the derived value of r, we plot in figure (2) the two constraints defined by equation 
(6) in the r] - Vf plane. The area between the curves defines the allowed region for the 
parameters. It can be seen that for any allowed value of the pulsed fraction a fraction of 
shots between 0.8 — 100% must be modulated at the spin period (coupled). Moreover, even 
in the most favourable case of 100% localized shots, the modulation has to be > 17% in 
order to account for the broadening detected. It can be noted that, in our discussion, we 
have assumed an identical amplitude for the localized and diffuse shots, that is, we have 
assumed that the blobs originating in the disk survive intact to the polar caps. It might 
be argued that the magnetic field just scoops off a fraction of each blob. We examined, for 
example, a case in which just 10% of each blob is threaded to the neutron star surface. In 
this scenario we obtain that nothing much changes in the favourable case of total coupling 
(100%), whilst for the maximum value of the pulsed fraction (100%), a minimum fraction 
of shots greater than 45% must be modulated at the spin period. That is, if just a part of 
each blob is scooped off, then more blobs must be funneled to observe the same effect. 

A different constraint can be derived comparing the "observed" pulsed fraction value 
"P/obs with its theoretical value Vf as derived from eq. 6 above. To compute the value 
of the "observed" pulsed fraction we have folded the data corrected for the orbital motion 
modulo the spin period and obtained a pulse profile. Adopting the background value for 
this observation of 130 c/s (Wijnands & van der Klis 1998a ) we found Vf bs ~ 12%. 
We can note from fig. (2) that the pulsed fraction 0.17 < Vf < 1 and, recollecting the 
definition of the fraction of diffuse component £ = 1 — Vt ^ s we can pose a constraint on £, 
1 _ El^i < £ < i _ gtob. and go 0.294 < f < 0.88. We can now plot in fig (3) the diffuse 

-L ' / min 

fraction £, versus the fraction of coupled shots r\. 

4. Summary and Conclusions 

We have shown the presence of a broadening at the base of the harmonic line. This can 
be interpreted as the result of a coupling between a fraction (between 0.8% and 100%) of 
the aperiodic and periodic variability in SAX J1808. 4-3658. Moreover we have demonstrated 
that a fraction between ~ 30% and ~ 90% of the total emission arises from a component 
that is not modulated by the lighthouse effect (diffuse component). 
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Some more stringent constraints can be derived if we assume that the diffuse component 
originates from an accretion disk and the localized component derives from accretion onto 
the neutron star surface. As the spin period of SAX J1808. 4-3658 is known, we can compute 
the corotation radius for this system i. e. the radius at which the speed of a particle rigidly 
rotating with the neutron star equals the local Keplerian value, R co = 1.5 x 10 m 

l/3p2/3 

cm, where m is the NS mass in units of solar masses and P_ 3 is the NS spin in milliseconds. 
In our case we get R co = 30.8 km assuming a mass of 1.4M . For the accretion not to 
be centrifugally inhibited the inner rim of the disc must be located inside the corotation 
radius. Therefore we have i?disc < Rco- In the standard scenario of disk accretion onto 
magnetized neutron stars, the accreting matter forms a disk whose inner radius is truncated 
at the magnetosphere by the interaction of the accretion flow with the magnetic field of the 
NS. In this case the magnetospheric radius r m is a fraction <fi < 1 (an expression for can 
be found in Burderi et al. 1998 1 ; for L ~ 10 36 ergs/s we get <fi ~ 0.3) of the Alfven radius 
i?A defined as the radius at which the energy density of the (assumed dipolar) NS magnetic 
field equals the kinetic energy density of the spherically accreting (free falling) matter: 

R A = 2.23 x l^R-^m^^e 2 ' 7 ^' 7 cm (8) 

(see, e.g., Hayakawa 1985), where Rq is the NS radius, R^s, m units of 10 6 cm, m is the 
NS mass in solar masses, /Z26 is the NS magnetic moment in units of 10 26 G cm 3 , e is the 
ratio between the specific luminosity and the specific binding energy (L = e x GMM/R NS , 
G is the gravitational constant, M is the NS mass, M is the accretion rate), and L 37 is the 
accretion luminosity in units of 10 37 erg/s, respectively. In this context we have -Rdisc = 
r m . Since the accretion disc is virialized, its luminosity (averaged over the shot process) is 
(7 DF ) = 0.5 eGMM/r m . The total luminosity is (I DF ) + (I LC ) = eGMM/R m . Therefore 
we have £ = 0.5 Ri<ss/ r m- A lower limit for the inner rim of the accretion disk is r m > i?Ns 
and therefore we get the constraint £ < 0.5. This constraint is indicated as a dashed 
horizontal line in fig. 3. This immediately gives a lower limit on the fraction of coupled 
shots rj > 42% indicated ad the dashed vertical line in figures 2 and 3. This in turn limits 
the value of the pulsed fraction below 0.37 as indicated in fig. 2 by the dashed horizontal 
line. An upper limit can be derived from the constraint derived in the previous paragraph 
£ > 0.294 that gives r m < 1.7 -Rns- Adopting as a typical NS radius -Rns ~ 10 km, we 
get r m < 17 km, that is consistent with the requirement that accretion is not centrifugally 
inhibited. Moreover adopting <fi = 0.3 we can derive an upper limit for the magnetic field 
of the NS. Adopting the explicit espression for Alfven radius in r m < 1.7 -Rns we obtain 
1^26 < 5.1 4>q 7 / 4 R^mr x l^e~ x l 2 l}J^, where phi , 3 is the factor in units of 0.3. Taking 



V = 0.21a 4 / 15 no / 6 2 1 7 5 TO^ 142 / 945 [(L 3 7/e) 8/7 i?6 /7 At26 7 ] 4/135 > see the text for the definition of the symbols. 
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rn = 1.4, e = 1 and L 37 = 0.1 that are appopriate for our source, we get /Z26 < 1-5 R% ■ This 
is well within the range (1 — 5) x 10 26 G cm 3 that are the lower and upper limits derived by 
Chakrabarty (199) and di Salvo & Burderi (2002) respectively. 

This scenario implies that the aperiodic variability is generated in a region affected by 
the lighthouse modulation, close to the neutron star surface. Our result is in line with the 
coupling of the red noise discovered analyzing the power spectra of some high mass X-ray 
binaries (B97, LS97). Our present result constrains some of the models proposed for the 
origin of the aperiodic variability. Indeed a scenario like the one proposed by Aly & Kuijpers 
(1990) in which the aperiodic variability is generated by a shot noise process associated 
with magnetospheric flares, caused by magnetic reconnection all around the neutron star 
magnetosphere, cannot work for the origin of the aperiodic variability in SAX J1808. 4-3658 
in case the red noise component were coupled with the periodic modulation. On the other 
hand, inhomogeneities could arise in the accretion flow during the plasma penetration in 
the magnetosphere (e.g. via Kelvin- Helmholtz instabilities, see B97 for a more extensive 
discussion). The subsequent impact of this inhomogenous accretion flow funneled by the 
magnetic field lines onto the magnetic caps appears to be a viable mechanism for the origin 
of the red noise. 

This work was partially supported by a grant from the Italian Ministry of University 
and Research (Cofin 2001021 123_002). We warmly thank the referee for very useful and 
stimulating criticism. 
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Fig. 1. — The power spectrum of SAX J1808. 4-3658. Upper panel shows the measured 
broad band power spectrum and the simultaneous modelling of the red noise and harmonic 
components (solid line). The lower panel shows an expanded view of the region around the 
first fundamental harmonic: data have been binned by a factor 35 to increase the signal to 
noise ratio. Both the fit with and without coupling with the red noise have been plotted. 
The lower line is without coupling. (Note that in the case of no coupling the residual 
broadening at the base of the harmonic line is due to the convolution of the delta function 
representing periodic modulation with the Fourier transform of the box function of duration 
T that multiplies the signal because of the finite duration of the observation). 



0.01 0.1 

Fraction of coupled shots 



1 



Fig. 2. — The Vf vs r/ relation. The dashed vertical line indicates a lower limit on the fraction 
of coupled shots i] > 0.42 derived from the constraints of Fig. 3. This, in turn, poses an 
upper limit on the pulsed fraction Vf < 0.37. 
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Fig. 3. — The £ vs rj relation. The horizontal line indicates an upper limit for the fraction 
of diffuse component £ < 50%. This, in turn, poses a lower limit on the fraction of coupled 
shots i] > 0.42, as shown by the dashed vertical line. 
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Table 1: Observations log 
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21:56 


n 


30411- 


-01- 


-02-00 


13- 


A 

-4- 


•1998 


01:23 


13-4-1998 


02:22 


n 


30411- 


-01- 


-03-00 


16- 


-4- 


-1998 


17:02 


16-4-1998 


22:59 


y 


30411- 


-01- 


-04-00 


17- 


-4- 


-1998 


01:11 


17-4-1998 


04:40 


y 


Qri/i i 1 
oU411- 


ni 
-Ul- 


-Uo-UU 


1 Q 


-4- 


1 nnc 


no. £ a 


iq /i i one 


nn. 1 1 
Uy.l 1 


y 


30411- 


-01- 


-06-01 


18- 


-4- 


■1998 


12:27 


18-4-1998 


13:32 


y 


30411- 


-01- 


-06-00 


18- 


-4- 


■1998 


14:04 


19-4-1998 


01:02 


y 


30411- 


-01- 


-07-00 


20- 


-4- 


■1998 


20:52 


20-4-1998 


23:07 


y 


30411- 


-01- 


-08-00 


23- 


-4- 


■1998 


15:47 


23-4-1998 


23:19 


n 


30411- 


-01- 


-09-01 


24- 


-4- 


■1998 


15:51 


24-4-1998 


17:14 


n 


30411- 


-01- 


-09-02 


24- 


-4- 


■1998 


17:28 


24-4-1998 


23:21 


n 


30411- 


-01- 


-09-03 


25- 


-4- 


■1998 


14:17 


25-4-1998 


15:28 


n 


30411- 


-01- 


-09-04 


25- 


-4- 


■1998 


15:49 


25-4-1998 


21:40 


n 


30411- 


-01- 


-09-00 


26- 


-4- 


■1998 


15:55 


26-4-1998 


23:26 


n 


30411- 


-01- 


-10-02 


27- 


-4- 


■1998 


14:19 


27-4-1998 


15:41 


n 


30411- 


-01- 


-10-01 


27- 


-4- 


■1998 


15:57 


27-4-1998 


19:29 


n 


30411- 


-01- 


-10-00 


29- 


-4- 


■1998 


14:08 


29-4-1998 


18:56 


n 


30411- 


-01- 


-11-02 


2- 


-5- 


■1998 


01:14 


2-5-1998 


02:57 


n 


30411- 


-01- 


-10-03 


2- 


-5- 


■1998 


03:37 


2-5-1998 


04:38 


n 


30411- 


-01- 


-11-00 


3- 


-5- 


■1998 


18:59 


3-5-1998 


21:00 


n 


30411- 


-01- 


-11-01 


6- 


-5- 


•1998 


10:52 


6-5-1998 


14:00 


n 
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Table 2: Results of the fit 





Parameter 


Value 


Error 

(90% confidence level) 


WN 




1.994 


±0.001 


RN 


A 1 


1.7x 


±6 x lO" 1 




Pi 


1.1 x 10 1 


±0.6 x 10 1 




P2 


2.6 


±1 




P3 


1.6 x 1CT 1 


±1 x 10- 2 




Tl 


3.45 x 1CT 1 


±| x 10~ 3 




OL\ 


3.3 


± 3 X 10- 1 




r 2 


6.5 x icr 2 


±0.5 >< ID" 2 




«2 


2.0 


±0.2 




T3 


3.9 x 10~ 3 


±1. x 10~ 3 




«3 


1.5 


±0.1 




A 2 


7.3 x 10 1 


±3 




A 3 


3.4 x lO- 3 


±1 x 10~ 3 


Coupling 


r 


0.002 


±0.001 



